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ABSTRACT 



This thesis gives a number of approximations and bounds 
for the renewal function in an ordinary renewal process. 

Each approximation and bound is calculated for the uniform, 
gamma and hyperexponential distributions and compared with 
the renewal function for these cases. They are also calcu- 
lated for the log-normal distribution and compared with re- 
sults of the simulation of the renewal function. Results 
are tabulated and shown graphically. 
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I. INTRODUCTION 



Let {X^,X 2 ,...} be a sequence of independent identi- 
cally distributed non-negative random variables. Let the 
distribution function of X., i = 1, 2,... be F(x) . Let 



S « X. + X ? + 
n 1 Z 



+ X n , n > 1, 



the n^h partial sum, and S q = 0 . Define the random vari- 
able 

N(t) = Sup { n | S n < t }, 

which is the number of renewals occurring before and in- 
cluding t. The renewal function, M(t) , is defined as 



M(t) = E(N (t) ) , t > 0. 



( 1 ) 



Let F n (t) be the n-fold convolution of F with itself. The 
relationship between M(t) and F(t) is given by the follow- 
ing, 

oo 

M(t) = Y, F n (t) ’ 

n=l 

which leads to the well-known renewal equation, 

t 

M(t) = F (t) + J M (t - x) dF(x) . (2) 

o 

The purpose of this thesis is to investigate approxima 
tions and bounds for the renewal function, since Equations 
(1) and (2) can only be solved for special cases of the dis 
tribution function F. 



In Section II formulae are collected for various ap- 
proximations and bounds that have appeared in the litera- 
ture (Barlow § Proschan (1965) , Bartholomew (1963) , and 
Butterworth and Marshall (1971)). These formulae are 
applied to the uniform, gamma, hyperexponential and log- 
normal distributions. In Section III the results of com- 
putations are presented in both tabular and graphical 
forms for each approximation applied to each distribution. 



II. APPROXIMATION FORMULAE 



Throughout this thesis the following notation will be 

used : 

E[Xi] = \ , 

x 2 e [X? ] 

k - — 2-i- - 1. 

and t 

F e ft) = X J (1 - F (u) ) du , t > 0. 

o 

With this notation F g is an equilibrium excess distribution 
of a renewal process, and from Smith (1958), if F is non- 
lattice , 



V? f 4- N . - \ f * 1 " • r 1 ^ 

V. L ) /v «■* • . 



r ^ 

V ~ ,» 



This equation, together with (2) , plays an important part 
in the development of approximations to M(t) . 

1 . A Simple Lower Bound 

Butterworth and Marshall (1971) show that for any 
renewal process 

n n 

♦ Y, F k M - Y V F k-lM - F n (t >- 

k=l k=l 



n > 0, (4) 

where F 0 (t) = 1, and when n=0 both summations are assumed to 
be empty. The right-hand side of (4) gives an increasing 
sequence of lower bounds which converge monotonically on 
M(t) . The case n=l will be used in this thesis and we define 



A(t) = At - F e (t) , t > 0. (5) 



R 



2 . An Improved Lower Bound 



The following general lower bound was found by 
Barlow and Proschan (1965) . Define 

B(t) - p^tX ■ !• w 

Then B(t) < M(t) . Since B(t) = (A(t) /F (t) ) , it is clear 
that B is an improvement over A when F (t) < 1. 

3 . Upper Bounds When F is NBUE 

A distribution of a non-negative random variable is 
called NBUE (new better than used in expectation) if 



I 



MW] 

l-F(t) 



dx < I , 



( 7 ) 



for all t where F(t) < 1. Equivalently F is NBUE iff F e 
> F(t) for all t. For details concerning this and other 
classes of distributions frequently used in reliability 
theory see Marshall and Proschan (1970) . Butterworth and 
Marshall (1971) show that if F is NBUE, 



C+- '» 
K'-J 



n n 

M(t) < Xt + Y, F k (t) " Z F e* F k-l Ct) 7 n > °> C8) 

k=l k=l 

where again F (t) = 1, and when n=0 , both summations are 
taken to be empty. The right -hand-side of (6) gives a de- 
creasing sequence of upper bounds which converge monotonic- 
ally to M(t) . 1 The case n=2 will be used in this thesis. 



i 

If the inequality is reversed in (7), F is called N1VUE 
(new worsa than used in expectation). In this case (8) give 
an improved sequence of lower bourds over the sequence given 
by (4) . 



We define 



C(t) = At + F (t) - F e (t) (9) 

= A(t) + F (t) . 

4 . Upper Bound When F is IFR 

A distribution is said to have increasing failure 
rate (is IFR) if Log(l-F(t)) is concave in t. Under this as- 
sumption Barlow 5 Proschan (1965) showed that if we define 



Dft'i - ^tF(t) 
D(t) ’ 



( 10 ) 



then D(t) > M(t) . 

5 . Approximations for M(t) 

Bartholomew (1963) derived the following approxima- 
tion for the renewal density m(t) - (dM(t)/dt) when F has a 
density f; 



{ . n ~ + . A F(t) 2 

m(t) ~ f(t) + F . 

G 

The following integral of this approximation is used as an 
approximation for M(t) , 



t 

E(t) = F(t) + A J du. (11) 

o e 

Another form of the renewal function can be derived from (1) 

or (2) . One can show that 

t 

M (t) = At + C [1 - F 0 (t-u)J dM(u) - F e (t). (12) 
o 

If we approximate dM(u) on the right hand side by Adu we get 
another approximation for M. Thus define, 



( 13 ) 



G (t) = At + A J [1 - F e (u)] du - F e (t). 
o 

Recall that A(t) and B(t) gave simple lower and upper bounds 
on M for F restricted to NBUE . However, (A(t) -At)-*--l and 
(B(t)-At )->0 as t-^ 0 . From ( 3 ) for F non-lattice it is known 
that (M(t) -At)-*k, so define 



H (t ) = At - F e (t) + ( 1 +k) F ( t) , 



( 14 ) 



the convex combination of A and B with the correct asymptotic 
behavior . 

6 . Distribution Tested 

The following distributions with specified param- 
eters are used to compare the bounds and approximations for 

M(t). 

a. Uniform Distribution (A=l) 



Then 



F(t) = \ 
= 1 



F e (t) = t - J- 
= 1 



These give: 



A(t) = ^ 



= t - 1 



if 0 < t < 2 
if 2 < t . 

if 0 < t < 2 
if 2 < t. 

if 0 < t < 2 
if 0 < t, 



B (t) = 4^ 

= t - 1 



if 0 < t < 2 



if 2 < t, 
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C(t) = 


ii + 1 

4 2 


if 0 < t < 2 


= 


t 


if 2 < t. 


D C t ) = 


2t 

4-t 


if 0 < t < 2 


= 


t 


if 2 < t, 


Eft) - 


4 log ( 4 ? t ) - \ 


if 0 < t < 2 


= 


t - 41og ( \ ) - 3 


if 2 < t , 


G(t) = 




if 0 < t < 2 


= 




if 2 < t, 


H (t) = 


ii + 1 

4 3 


if 0 < t < 2 


_ 


t - \ 


if 2 < t. 




3 




For this case. 






M (t) = 


V 1 r t i n 2' 1 ' 

L nT ("T’ e 


-1 if (j-l)(2) < t < ( j ) (2) 




_ n=0 


J 1 y 2 y • • • • 


b. Gamma Distribution (A=l) 


F(t) = 


1 - (1 + 2t) e " 2t , 




- 


1 - (1 + t) e" 2t . 




These give : 






A(t) = 


t + (1 + t) e' 2t - 1 , 


B(t) = 


^ -1 

_ 9 + - 1 - 9 




1 - (1 + t) e 





t 



C(t) = t - t e' 2t , 



D (t) = 



t-te 2t -2t 2 e 

" I Tt 7 r 2T 

1 - e - t e 



2t 



E(t) = 1 - e' 2t - 2 t e' 2t + 



/‘H 



- 2u - - 2u^ 2 

e -2u e ) , 

e" 2u - u a’ 2 " ’ 



G (t) = t + ( \ t + ) e‘ 2t - 

H(t) = t - \ ( t - \ ) e' 2t - j. 



For this case 



X(ir .s . . 1 -4t 1 

M(t) = t + j e - 2 J-. 

c. Hyperexponential Distribution (A=l) 



F(t) -4(1- e' 2t ) + t ( 1 - e b ) , 



F e (t) = | ( 1 - e “ ) + £ ( 1 - e " ). 



-2t 



5 

8 






These give: 



A (t) = t ♦ |e' 2t + |e'5 * 



- 1 , 



B(t) = 



8t 



- 2 t 

3(1 - e ” 2t ) + 5(1 - e 5 ) 



C(t) = t + | ( 1 - e' 2t ) - | ( 1 - e 5 t ) , 



-2 t 

2t [3 (1 - e _2t ) + ( 1 - e 5 )] 

D(t) - 2 » 

t 

3 ( 1 - e‘ 2t ) + 5 ( 1 - e 5 ) 



E(t) = | ( 1 - e " 2t )*J(l-e st ) 



-2u ‘t u 

+ , 1 3(1 - e ) *J 1 ^ e ,5 du _ 

9 “F U 

0 6(1 - e~ zu ) + 10(1 - e 5 ) 



t 

■i 



G(t) = t - T# ( 1 - e' 2t ) * i| ( 1 - e 5 ), 



H(t) = t 



16 

15_ 

16 



C 1 - e" 2t ) 



16 

3 



-ft 

TZ ( 1 - e ) 



For this case 



Then 



3 4* 

M(t) = t + J ( 1 - e b ) 



d. Log-normal Distribution (A=/e — ) 
Let 



<Kt) = 



t XJ 
2 



— f * 

/2¥ J 



dx . 



F (t) = $ (log t) 



F P (t) = te (1 - 3>(log t)) + $ ((log t) - 1) 



1 

2 



These give : 

A(t ) = e t <Hlog t) - $ ( (log t) - 1) , 

B (t ) = 



1 

e' 2 t 



e 2 t( 1 - 



$(log t)) +$ ((log t) - 1) 



C(t) 

D(t) 



E(t) 



_1 

(e 2 t + 1) <Hlog t) - $((log t) - 1) , 

_1 

e 2 t $(log t) 

_1 

e 1 t Cl -$Clog t)) + $((log t) - 1) 

.1 r 

*(log t) + e 2 I — T Cniog u)) 2 du 

o e 2 u(l-$(log u))+$((log u) -1) 



_1 

G(t) = e 2 t - | Cl - $C Clog t) - 2)) 

_1 

+ Ce 2 t + 1) Cl " <K(log t) - 1)) 

_1 _1 

- \ e 2 t Ce 2 t + 2) Cl - <Hlog t)) + | - 1, 
_ 1 _ 

HCt) = C e 2 t + \ e) $Clog t) - $(Clog t) - 1) . 
For this case the exact form of M(t) cannot be found. 



III. NUMERICAL RESULTS 



Numerical calculations are presented for each bound and 
approximation and all four distributions in Section II. Re- 
sults are tabulated below and are graphed and comparisons 
are made with M(t) . For the log-normal case M(t) was simu- 
lated for comparison with approximations. 
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FIGURE 1 : M(t) and E(t) for the uniform distribution. 



G(t) 
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FIGURE 2 : M(t) and G(t) for the uniform distribution 
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FIGURE 3 : M(t) and H(t) for the uniform distribution 
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TABLE II. M(t) and approximations for the gamma distribution with mean 1 and 
variance 1/2. 
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FIGURE 4 : M(t) and E(t) for the gamma distribution. 
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TABLE III. M(t) and approximations for the hyperexponential distribution with 
mean 1 and variance 5/2. 
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TABLE IV. M(t) and approximations for tlie log-normal distribution with 
and variance (e 2 - e) . 
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FIGURE 11 : M(t) and G(t) for the log-normal aistributi 
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FIGURE 12 : M(t) and H(t) for the log-normal distribution. 
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